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Electronic spin systems with S > 1/2 provide an efficient method for DC vector magnetometry, 
since the conventional electron spin resonance spectra at a given magnetic field reflect not only the 
field strength but also orientation in the presence of strong spin-spin interactions. S=1 spins, e.g. 
the nitrogen-vacancy centers in diamond, have been intensively investigated for such a purpose. In 
this report, we compare S=1 and S=3/2 spins, and discuss how one can apply general principles 
for the use of high spin systems as a vector magnetometer to the S=3/2 spin systems. We find 
analytical solutions which allow a reconstruction of the magnetic field strength and polar angle 
using the observed resonance transitions if an uniaxial symmetry exists for the spin-spin interaction 
as in S=1 systems. We also find that an ambiguity of determining the field parameters may arise 
due to the unique properties of S=3/2 systems, and present solutions for it utilizing additional 
transitions in the low-field region. The electronic spins of the silicon vacancy in silicon carbide will 
be introduced as a model for the S=3/2 DC vector magnetometer and the practical usage of it, 
including the magic-angle spinning type method, will be presented too. 


I. INTRODUCTION 

Electronic spins in highly localized defects, such as 
the nitrogen-vacancy (NV) centers in diamond^’^ and va¬ 
cancy related defects in silicon carbide may 

experience a strong spin-spin interaction, e.g., a dipole- 
dipole interaction, which results in the so-called zero-field 
splitting (ZFS), partially (or completely) lifting degener¬ 
acy of energy eigenstates at zero magnetic field^^’^®. If 
this interaction is strong enough, the eigenvalues of the 
spin Hamiltonian show a strong dependence on the ori¬ 
entation of the applied magnetic field. Such dependence 
causes a non linear shift of resonance transitions in elec¬ 
tron spin resonance (ESR) spectra. Thus, the informa¬ 
tion about the applied external magnetic field can be 
extracted from ESR spectra provided the ZFS is known. 

One well-known example is the NV center in diamond. 
Its application to DC field vector magnetometry has been 
reported and well understood in the field strength from 
sub-/xT to a few tenth T^^d 8 . 20 , 2 i^ 
a triplet ground state of S = 1 and when shifts of the 
ESR transition at a given DC magnetic field are directly 
monitored in the frequency domain, typically a ~ 0.1 mT 
minimum detectable magnetic field is achieved^®. This 
resolution is limited by the ESR linewidth which can be 
broadened by strong RE fields thus lowering the resolu¬ 
tion. Lower RE power can be used to avoid power broad¬ 
ening, but the decreased signal strength requires a very 
long accumulation time. If time-domain experiments, 
e.g. a Ramsey fringe experiment, in which the magnetic 
field strength is imprinted in the phase of the superpo- 
sitioned state, is conducted, a large signal strength can 
be maintained without power broadening, thus a sensi¬ 
tivity up to ~ 0.4 /xT/VHz, limited by the TJ of ~ 1 /xs, 
can be realized using a single NV center^’^°. Further en¬ 
hancement (below 1 nT/-\/Hz) is possible by using the 
NV center ensemble combined with the lock-in detec¬ 
tion^^. When the NV center is used for AC magnetic 
held sensing, spin echo type measurements can be used in 


which the long coherence time allows high sensitivity up 
to ~ 1 nT/VHz using a single NV^^ and ~ 0.9 pT/-\/Hz 
using an NV ensemble^^. 

Higher spin systems (S > 1) can also be used as a vec¬ 
tor magnetometer in a similar way. For example, the 
silicon vacancy (Vsi) in silicon carbide (SiC) is known 
to possess a quartet manifold of S = 3/2 in its electronic 
ground state^^“^®. Because its ESR signal can be de¬ 
tected at ambient condition^’^^’^^“^® even from a single 
defect^ and the ZFS is in a range around a few mT de¬ 
pending on the polytype of SiC^’^, its application as a DC 
magnetometer has been suggested^®’^®. Note that Simin 
et al, have recently shown an experimental application 
of Vsi in SiC as a sub-mT DC magnetometer based on 
approximated solutions for the spin Hamiltonian at weak 
magnetic helds^®. 

In a spin system with a spin quantum number S, the 
strength and orientation of the applied magnetic held 
vector Bq determines the Zeeman splitting, thus one 
should experimentally obtain the Zeeman splitting to get 
information about Bq. For S = 1/2, the Zeeman splitting 
is calculated from an observed single resonant transition 
energy hi/ = g^^Bo where g is the Lande g-factor, fig 
is the Bohr magneton, and h is Planck’s constant. The 
information for the orientation can be extracted only if 
g is anisotropic. In high spin systems, the orientation 
related terms remain in the eigenvalue equation which 
result in the orientation dependent shift of ESR spec¬ 
tra, which cannot be explained by g/isHo- K is, there¬ 
fore, mandatory to reconstruct the energy eigenstates us¬ 
ing the observed resonant energies. Because there exist 
2S-I-1 eigenstates, 2S resonant transition energies should 
be experimentally determined. For example, when the 
applied magnetic field strength is much larger than the 
ZFS, i.e., (//xb|Bo| ZFS, at least two transition en¬ 
ergies should be known for S = 1 while three values are 
necessary for S = 3/2 as explained in Fig.I. In this high 
field range, the two transitions for S = I and two out 
of the three transitions for S = 3/2 cross each other as 
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FIG. 1. (Color online) Energy eigenvalues and ESR transition 
energies of high spin systems as a function of orientation at 
a high magnetic field gpsBo = 100 x ZFS. (a) and (b) are 
the eigenvalues of the spin triplet (S=l) and quartet (S=3/2) 
state calculated numerically from Eq.(l), respectively. For 
numerical calculation, D,E > 0, and an uniaxial symmetry 
E D are assumed thus only the polar angle, 6 dependence 
is shown. Eigenstates are labeled in ascending order of corre¬ 
sponding energy eigenvalues. The most dominant transitions 
are indicated by the solid arrows with labels fij indicating a 
transition between |i) and |j). They are shown in (c) and (d) 
for the spin triplet and quartet states, respectively. All the 
energy values are normalized by ZFS. 


shown in Figs.1(c) and 1(d). This leads to ambiguity 
in determining which observed ESR peak corresponds to 
which transition energy experimentally. In this report, 
we discuss how this ambiguity can be removed and thus 
show how to use the S = 3/2 system for vector magne- 
tometry. For this, we will provide analytical solutions for 
the given Bq vector as a function of the resonant transi¬ 
tion energies. We will present Vsi spins in SiC as a model 
system and also a novel magnetometry scheme using the 
magic angle. The discussion in this report is also appli¬ 
cable to other S = 3/2 systems which have been found 
in fullerene^®^^"^, organic molecules^®^^’^, Ni impurities in 
diamond^®, and calcium oxide crystals®®. 


II. VECTOR MAGNETOMETRY BASED ON 
S=3/2 SPINS 

In order to derive formulas for Bq = |Bo| and its ori¬ 
entation expressed by only three transition energies for 
S = 3/2, we will first construct the electronic spin Hamil¬ 
tonian consisting of the ZFS and Zeeman term. The ZFS 
in high spin systems can be described by the dipole-dipole 
interaction term in the spin Hamiltonian, S • D • S where 
D is the dipole-dipole coupling tensor. For simplicity, 
we assume an isotropic Lande g-factor. Therefore in the 
principal axis system of D, in which the z-axis is set to 
the symmetry axis, the electronic spin Hamiltonian at 
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FIG. 2. (Golor online) Bq dependence of the energy eigen¬ 
values of a spin quartet state at various magnetic field ori¬ 
entations. Eigenstates are labeled in ascending order of the 
corresponding energy eigenvalues. As Rq oo at 0 = 0°, 
|4) ^ |ms = +3/2), |3) ^ |ms = +1/2), |2) ^ |ms = -1/2), 
and |1) —>■ |ms = —3/2). 

Bo, is 

H = g^iBBQ-S+D{Sl-S{S+l)/3}+E{Sl+Sl)/2, (1) 

where E and D are the ZFS parameters, assumed to be 
positive, and E <?; H if an uniaxial symmetry exists. For 
S=3/2, the eigenvalue equation from Eq.(l) is, in the 
polar coordinate system, 


- (2H® + 6E;2 + - 2l3l{D{3 cos^ 6 » - 1) + dE sin^ 9 cos 2(/}A 

+ - T'®/3^(3cos2 9-1)+ 3E^{3E^ + 2D^) 

16 2 

+EI3q {6D sin^ 9 cos 2(j)+^E cos 29) = 0, (2) 


where /?o = gg-sBo- The numerically calculated eigen¬ 
values at various orientations are shown in Fig.2. When 


Bo is either parallel or perpendicular to the symmetry 
axis, the closed form solutions for each eigenvalue can be 
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found as^® 


A = ^/3o ± \/{D + /3o)^ + or - ^j3o ± ^(£> - /3o)^ + for Bo||z - axis, 

A = ^/3o ± y ^/32 + i ^2 + 3^2 _ or 

-^/?o ± \JPo + + (D - 3£^)/3o for Bq || x - axis, (3) 


which gives the eigenvalues at zero magnetic field, 
Xbo=o = ±ZFS/2 where ZFS = 2VD^ + 3E^. The eigen¬ 
value equation for the general case can be expressed as, 

2S-I-1 

CnX^ = 0 . ( 4 ) 

n—0 


By plugging each eigenvalue Xi into Eq.(4), one can ob¬ 
tain 2S-I-1 equations. The basic idea in order to find 
formulas for the Bq vector expressed by the observed res¬ 
onant energies, is to remove all Xi terms using the tran¬ 
sition energy = Xi — Xi-i. Note that the energy 

eigenstates are not necessarily sorted with respect to the 
corresponding energy values. In other word, the indices 
can be randomly assigned to the states. Here, however, 
we keep the relation, Xi > Ai_i, for convenience. We 
follow this approach which has been frequently used for 
S=1 systems^^’^®. First, (2S-1) sets of three simultane¬ 
ous equations are obtained by plugging Xi -I- fi+i^i, Xi, 
and Xi — fi^i-i (i=2,3,..2S) into Eq.(4). In each set, cal¬ 
culating 


2S-I-1 

E 


n—0 


Cn{{Xi + 


C 2 S+I 


m 


= 0 and 


2S-I-1 

E 


n—0 


C„{(A, 


1^25-1-1 



( 5 ) 


results in two new simultaneous equations 


2S 2S 

5 ]C',„Ar=0and ^ 


^ C”nK = 0. 


( 6 ) 


n—0 


n—0 


Again, by subtracting one from each other divided by 
the coefficient of the highest order term of each equation, 
respectively, as below 


2S 

E 


/^/ \ n \ V 


r" 

n—0 ^,25 


c: 


= 0 , 


( 7 ) 


i,2S 


we obtain a new equation for the eigenvalue of the energy 
eigenstate |t) in which the highest power is (2S-I), 


2S-1 

E = 0 , (8) 

n—0 

where i=2,3,..2S. By repeating this procedure until only 
one linear equation for a single eigenvalue remains, one 
can find a formula for an eigenvalue expressed in terms 
of resonant energies, which allows us to find expressions 
for all other eigenvalues again using We use this 

procedure to find solutions for S=3/2. 

Following the procedure explained above, we obtain 
two equations for two energy eigenvalues, A 2 and A 3 , 
expressed by only / 2 ,i and /a, 2 , and /a ,2 and 74 , 3 , re¬ 
spectively, and Bq, 9, and (j) which are present in both 
equations. Then using /a 2 = A 3 — A 2 once again, we ob¬ 
tain formulas for each eigenvalue expressed by only the 
resonant energies as below. 


^1 = “4/2.1 “ 2/3.2 “ 4/4.3, 

-*^2 = 4/2.1 “ 2/3.2 “ 4/4.3, 

-^3 = 4/2.1 + 2/3.2 “ 4/4,3, 

A 4 = ^/2,i + 2 / 3.2 + 4 / 4 , 3 - (9) 


And by plugging one of these, e.g. A 2 , back into one 
of the equations found in the preceding steps, we finally 
obtain formulas for /3 q, and a new quantity related to 0 
and (j), rj = E{2cos^(j)sm‘^6 -|- cos^9) + Dcos^ 9, 


,2 (#/4,3 + /3.2 + #/2,l)" + (1 - 73)(/4.3 + /2.l)/3.2 - /4.3/2.I - ZFS^ 

0 5 

4[8(iA + 3E) + 5 (/ 4.3 - f2,iM + (/ 4.3 - / 2 .i)[ 4 ZFS^ - (A ,3 - / 2 . 1 )" - 4 / 3 %] 
^ 96/32 


( 10 ) 


( 11 ) 
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FIG. 3. Orientation dependence of ESR transitions at (a) 
high (jr/iflBo = 100 x ZFS) and (b) low (qubBo = ZFS/10) 
magnetic field strength. The dashed lines indicate the magic 
angle 6m.- The gray color scale depicts the normalized ESR 
transition probabilities for Bi|| x — axis 


Thus, in general, if D and E values are known and three 
resonant energies are observable, the applied magnetic 
field strength can be extracted from Eq.(lO). In addi¬ 
tion, if an uniaxial symmetry exists (£1 <C D), because of 
Tj D cos^ 9, the polar angle can also be extracted from 
Eq.(ll). 

However, when one wants to use these formulas to find 
the external magnetic field vector, one may face a prob¬ 
lem: how can one determine which observed resonant 
energy is from which transition? This is explained in 
Fig.l and Fig.3. At a given magnetic field vector (e.g., 
gfJ.BQ ^ ZFS), there will appear three strong transi¬ 
tions for S = 3/2 as in Fig.3(a). The resonant ener¬ 
gies are varying depending on the orientation for fixed 
Bq and some of them even cross each other, thus it is 
hard to determine fij explicitly only taking into account 
these three transitions. However, in certain spin systems 
bound to certain localized defects in solids and at a cer¬ 
tain magnetic field range, this ambiguity can be relaxed 
as will be shown in the following sections. We will fo¬ 
cus on two distinguishable cases for gg-sBo <C ZFS and 
ggsBo ^ ZFS. The case for ggsBo ~ ZFS, however, 
will not be discussed because complex spectra appear 
due to a strong interaction among each eigenstate^® via 
e.g. the level anti-crossing as in the NV centers^^, thus, 
high spin systems are not an appropriate sensor for this 
field range. In the following sections, we will consider 
the uniaxial symmetry for all the numerical simulations 
unless noticed. This will allow us to utilize Eq.(ll) to 
extract at least one orientation component, 9. 

A. High field, gfisBo S> ZFS 

Figure 3(a) depicts how each transition evolves at vary¬ 
ing orientation 9 at high static magnetic fields {ggsBo = 


100 X ZFS) together with ESR transition probabilities. 

is assumed to have only an x-axis component. At 
high fields, three transitions are visible. Without know¬ 
ing the information about Bq, it is not possible to as¬ 
sign them correctly because as will be seen later, there 
appears also three or even more than three resonances 
that originate from different transitions at low magnetic 
fields. However, one can guess Bq roughly from the ob¬ 
served resonant energies because they are approximately 
proportional to Bq at high field. For example, / 3^2 shows 
very weak orientation dependence thus can be used to 
estimate Bq. Once the Bq scale is roughly guessed, one 
can try to assign the observed resonant energies. Be¬ 
cause /a ,2 always stays in the middle of all resonance 
transitions, this can be explicitly determined. And be¬ 
cause Eq.(lO) is invariant under switching of and 
/ 2 . 1 , Bq can be unambiguously extracted. In contrast, 
Eq.(ll) will result in a systematic error if / 4_3 and f 2 ,i 
are not correctly assigned. Thus, a magnetometer based 
on a S' = 3/2 system can be used only to extract Bq at 
a high magnetic field. This is a disadvantage compared 
to S = 1 because a similar equation to Eq.(ll) is also 
invariant under switching of two observed resonant en¬ 
ergies'^. This problem, however, may be overcome by 
manipulating the ZFS^®. If the electric dipole moment is 
large enough and their contribution to the ZFS Hamilto¬ 
nian, the Stark-effect, is well-known, manipulation of the 
ZFS by either applying an electric field^^ or pressure"^^ 
along a favored direction can result in the shifts of each 
ESR transitions in different manners. Thus monitoring 
the additional shift upon the change in ZFS, may allow to 
determine all the necessary transitions and subsequently 
Eq.(ll) can be used without systematic errors. 


B. Low field, gfisBo ^ ZFS 


Figure 3(b) shows that one can see up to five tran¬ 
sitions at low field. At a small angle, there are three 
dominant transitions, / 2 p, and / 4 _ 2 , and two ad¬ 
ditional transitions, / 3 , 2 , and / 4 _i, which arise at large 
angle. Alternate forms of Eq.(lO) and (11) using only 
the most dominant transitions can be found for small 
angles (not shown). These forms, however, are not use¬ 
ful because / 4^2 and are changing their relative po¬ 
sitions at a larger angle. Instead, / 4 ^i and / 3^2 can be 
used since their relative positions are not changing and 
in certain S = 3/2 systems, e.g. Vsi in SiC (see section 
HI), these transitions show good intensities at every ori¬ 
entation^®. The other useful formulas can be found by 
plugging / 4,3 = / 4 q - / 3,2 - / 2 ,i into Eqs.(lO) and (11) 
as 
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q 2_ (y^favg + /2,l)^ — 2/4^i/3_2 + (1 ” ■\/3)/3, 2 / 2,1 ” (1 + ■\/3)/4l/21 — ZFS^ 

Pq— ^ 

^ + 3E) + 5Afout - Wf2M - 8/2.i(ZFS2 - fl, - fl^) 

yopp 

+A/o„i[4(ZFS2 - 3fl, - /3%) + A/o„i(6/2.i - Afout)]}, (13) 


where Afout = / 4.1 - /a ,2 and favg = (/s.i + / 4 . 2)/2 = 
(/a ,2 + /4,i)/ 2- These are useful, because if and / 3_2 
are observable together with / 4_2 and / 34 , one always can 
unambiguously determine the two outermost transitions, 
and / 4^2 and are not in use or necessary only for 
calculating favg- 

So far, the strategies to use S = 3/2 systems as a 
DC vector magnetometer have been discussed in both 
high and low magnetic field ranges. Though only Bq 
can be obtained at high fields using only a conventional 
experimental method, both Bq and polar angle can be 
determined at low field. However, in many high spin sys¬ 
tems bound to localized defects in solids, spin-dependent 
intersystem-crossing may induce a strong polarization 
into specific spin states as in Vsi in SiC^°’^®. Thus, some 
transitions may be hardly observable. In addition, be¬ 
cause only the polar angle can be obtained, it is still not 
possible to realize a genuine vector magnetometry. In 
the following sections, we will present S' = 3/2 spins of 
Vsi in SiC as a model system and discuss the practical 
usage of them and a possible way to use them as a vector 
magnetometer. 

III. SILICON VACANCY SPINS IN SILICON 
CARBIDE AS A DC VECTOR 
MAGNETOMETER 

We present Vsi in SiC as a model system to provide ex¬ 
planations about how the formulas, found in previous sec¬ 
tions, can be used to experimentally reconstruct the ap¬ 
plied external magnetic field vector. Because spin prop¬ 
erties are different depending on the polytype, here we 
discuss only a specific polytype, namely 4H-SiC. In ad¬ 
dition, because there exist two inequivalent lattice sites, 
there appear two different silicon vacancies with differ¬ 
ent ZFS, and we choose only one of them known as Tv 2 a 
center^’^. In the Tv 2 a center, it is known that there exist 
an uniaxial symmetry around the c-axis thus E D and 
ZFS/ft. ~ 2D/h ~ 70 MHz^’"*. It is also known that opti¬ 
cal polarization results in equal populations in two sub¬ 
states, I ms = ±1/2). This is responsible for the absence 


of a transition between |ms = ±1/2) and |ms = —1/2) 
while aother two transitions, between |ms = ±3/2) 
and |ms = ±1/2), and between |ms = —3/2) and 
I ms = —1/2) are observable in the ESR spectra of Tv 2 a 
for Bo||c-axis^’®’®“^^’^^“^®. In the Bq orientation depen¬ 
dence at low and high magnetic fields^’^®’^’^’^®, one of 
the allowed transitions, corresponding to the transition 
between |ms = ±1/2) and |ms = —1/2) for Bo||c-axis, 
has not been observed probably due to that this equal 
population is somehow maintained. This will prevent 
Eq.(lO), (11), (12), and (13) from being used because 
fz ^2 at high fields and /2,i at low fields will not be observ¬ 
able. This transition, however, can become visible once 
electron-electron double resonance (ELDOR) is applied. 
The population difference between |ms = ±1/2) states 
can be induced by applying e.g. a resonant tt pulse be¬ 
tween \ms = ±3/2) and \ms = ±1/2) (or \ms = —1/2) 
and \ms = —3/2)) states which enable detection of this 
missing transition^^. This will allow unambiguous deter¬ 
mination of one transition /2p at low fields or /3^2 at high 
fields experimentally. 

At high magnetic field (e.g., Bq ~ 300 niT) as in 
Fig.3(a), two outer transitions, /4_3 and /2,i have been 
observed experimentally at almost all orientations at 
both cryogenic and room temperature^® except the 
central peak. The central peak is observable by ELDOR 
experiments®®, thus Eq.(lO) can be used as explained in 
section (II A). However, the polar angle, 9, can be deter¬ 
mined from Eq.(II) only at small angles because of the 
ambiguity on determining 74^3 and /2,i at larger angles. 

ESR spectra of Tv 2 a centers at a low magnetic field 
(e.g. sub-mT) as in Fig.3(b) allow an unambiguous de¬ 
termination of both Bq and polar angle as long as the 
ELDOR can be used to determine f 2 ,i as explained in 
section (HB). However, one can consider another case 
in which either ELDOR experiments are not available 
or /2,i is hardly observed in the ELDOR spectrum. In 
such a case, if /4p and /3_2 are observable, using rela¬ 
tions -I- /4_2 = /a,2 ± /4,i) /4,a = /4,2 — fa, 2 , and 
/2,i = /a,i — fa, 2 from Eq.(9), we again obtain alterna¬ 
tive forms of Eqs.(IO) and (II) as 


,2 _ {Vifavg + /a, 2 )" - / 4 . 2 /a.l - 2(V3 + l)fa,2favg - ZFS® 
° 5 

[32{D + 3E) + 20A/„]H2 + A/,„(4ZFS® - A/® - 4/3®2) 

96B^ 


( 15 ) 
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where A fin = — f 3,1 ■ Note that / 3_2 appears in both 

formulas but because it is always the lowest energy tran¬ 
sition, this can be explicitly determined. Similarly, one 
can find additional alternatives using f^ i instead of / 3 ^ 2 - 
Therefore, even if ELDOR is not available, as long as ei¬ 
ther / 4 ,i or / 3^2 is observable together with / 4_2 and /s^, 
Bq can be extracted using Eq.(14) because it is invariant 
under switching /4 2 and /sp. This scheme is feasible 
since / 4 p and / 3 p are observable from Tv 2 a in SiC by 
cw methods with a decent signal strength at sub-mT as 
recently reported^®. Eq.(15), however, still cannot pro¬ 
vide an unambiguous way to determine the polar angle 
because of A/i„ which changes signs if / 4 p and / 3 p are 
not correctly determined. 

So far, the strategies to use Vsi in SiC as a vector mag¬ 
netometer has been discussed. While the magnetic field 
strength can be extracted in both high and low magnetic 
field range, the orientation can be extracted only if there 
exists an uniaxial symmetry at a low magnetic field, and 
the azimuthal angle cannot be determined in any case. 
Note that the S=1 system with the uniaxial symmetry 
also can provide only the polar angle. But in the case of 
the NV center in diamond, because the NV centers can 
be in four different orientations along the diamond bond 
axes, one can determine both the polar and azimuthal 
angles from the shift of transitions of the inequivalently 
oriented NV centers. Similarly, Vsi in inequivalent lat¬ 
tice sites, e.g. Tyia and Tv2a in 4H-SiC, and Tyiaj Ty2a 
and Ty 3 a in dH-SiC^ can also be utilized. However, ESR 
spectra of Tyia and Ty 3 a are hardly visible at room tem- 
perature^’®’^^. Thus, an alternate method relying only on 
the Tv 2 a center that can be used for any magnetic field 
strength at room temperature is necessary. In the next 
section, another method using a magic angle that allows 
for the use of S=3/2 as a vector magnetometer will be 
discussed. 


IV. VECTOR MAGNETOMETRY USING 
MAGIG ANGLE 


We start from the eigenvalue equation in Eq.(2). In 
this equation, one can find terms including (3cos^0 — 1), 
which becomes zero at the magic angle 9^ — 54.7°. 
Eq.(2) can be simplified for 9 = dm and E D as, 

A 4 _ (2^2 + + I/34 + ^ 0, (16) 

and the eigenvalues are simply 

X = ±yAD^ + 5/32 ± 4 ^ 3/32742 (17) 

as depicted in Fig.2(c). For high Bq, these can be again 
approximated as 

^ = or ± ^^0 (for g^BHo>ZFS). (18) 




frequency (MHz) 


FIG. 4. (Color online) Magic angle DC magnetometry based 
on a S=3/2 spin system, (a) describes the frames used for 
the numerical simulation. See text for details, (b) shows the 
numerical simulations of the cw ESR spectra when the S=3/2 
spin is misaligned hy 6 = Om relative to the rotational axis at 
various orientations of Bo vector in y' — z' plane at low field 
Bo= 6 G. For the simulation, the ZFS of Tv 2 a center in SiC, 
ZFS/h ~ 2D/h ~ 70 MHz is assumed, (c) shows the same 
simulation but for the case when 9 6m and 0 = 30 °. 


Thus, at 9m, we obtain |Ai| = IA 4 I and IA 2 I = IA 3 I, and 
can see the least number of transitions as seen in Fig.3. 
We can use this aspect to use S = 3/2 system as a vector 
magnetometer. If the spin sensor is being rotated around 
an axis and the orientation between the c-axis and the 
rotational axis is fixed to 9m, one expects to see the least 
number of transitions whose widths are the narrowest 
when the rotational axis is aligned to the applied external 
magnetic field. In contrast, when the rotational axis is 
misaligned, very broad ESR transitions appear due to 
orientation sweeping. 

Figure 4 describes such an experiment in which the SiC 
crystal is attached to a rotational axis forming 9m with 
respect to the c-axis of the crystal. The resonant RF can 
be applied using a miniature coil surrounding either the 
rotational axis or the SiC crystal, similar to what has 
been suggested for the quantum gyroscope based on the 
NV center^"^. For the detection, a small sized ESR cav¬ 
ity can be used for conventional ESR detection. Fiber 
coupling also can be considered for ODMR. If an electri¬ 
cally detected magnetic resonance is possible, which has 
recently been shown in high spin systems'^® and also in 
SiC^®, additional small circuits can be utilized. Bi field 
modulation can be used to enhance the signal to noise ra- 
tio^^ because the ELDOR experiment which requires Bi 
pulses is not necessary for this experiment. For the sim¬ 
ulation of ESR spectra in this experiment, a laboratory 
frame is assumed: the z' — axis is set to the rotational 
axis. The angle (0) between the c-axis of the crystal, 
rotating with a constant speed around the z' — axis, and 
the external magnetic field can be derived and replace 
9 in Eq.(2). For convenience, the RF field is assumed 
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to be in the x-axis of the rotating frame. By assuming 
a Lorentzian lineshape with 3 MHz FWHM, the numer¬ 
ically simulated ESR spectra at varying d' while 0' is 
fixed to 90° are simulated for a low field {Bq = 6 G) 
as shown in Fig.4(b). As expected, when 0' = 0 equiv¬ 
alently <3 = dm or Bq II z' — axis, the narrowest tran¬ 
sitions are found while seriously broadened peaks like a 
powder pattern appears when misaligned {9' ^ 0). Note 
that in order to present a general S = 3/2 case, / 2 ,i is 
assumed to be visible in cw ESR spectra at low field. 
Therefore by monitoring the linewidth of the observed 
transition spectra while moving the rotational axis, z'- 
axis, one can explicitly find the orientation of the external 
magnetic field. The field strength also can be extracted 
from the observed resonant energy of the strongest tran¬ 
sition using Eq.(17). Eor very small field gfisBo ^ ZFS, 
/ 3.1 = / 4,2 ~ gg-sBo + 'iggsBoD'^/3. We also can ob¬ 
serve spectra consisting of the narrow transitions even if 
9 ^ 9m 8i,s long as 9' = 0. However, because many transi¬ 
tions whose intensities are comparable to each other ap¬ 
pear as in Fig.4(c), it is more convenient to use the magic 
angle because the most dominant transition is easily dis¬ 
tinguishable. 

V. CONCLUSION 

We have shown that using the Vsi in SiC as a model 
system, S=3/2 electronic spins with the uniaxial symme¬ 
try can be used to find the strength and polar angle of 
the applied external magnetic field if at least three ESR 
transitions can be found experimentally and the ZFS pa¬ 
rameters are known. At a high Bq field {ggsBo ^ ZFS), 
Bq can be obtained from the observed ESR spectra but 
the polar angle cannot be determined due to the ambigu¬ 
ity of differentiating two outer transitions. In contrast, 
at low g/isHo (^ ZES), as long as one can explicitly 
identify at least three transitions including the allowed 
lowest energy transition, the external magnetic field vec¬ 
tor can be reconstructed. In the field strength compara¬ 
ble to the ZES, it is hard to find a useful scheme because 
very complex patterns appear due to mixing of some of 
the eigenstates. In the case of the NV centers in dia¬ 


mond (ZES/h=2.87 GHz), this missing range is around 
~ 100 niT. The Vsi in SiC can fill out this gap since its 
ZFS is quite small (ZFS/fi ~ 100 MHz) thus this mag¬ 
netic field range can be considered as a high field range 
in which the three necessary transitions are well observ- 
able^^’^®, and at least the field strength can be experi¬ 
mentally determined. When the Vsi in SiC is used to 
realize such schemes at sub-mT, if the lowest transition 
energy is observable by ELDOR, one can determine both 
Bq and 9 without ambiguity. Even if ELDOR is not avail¬ 
able, thanks to the additional transitions that appear at 
low fields, the field strength can be determined. 

The magic angle terms in the eigenvalue equation al¬ 
low for an alternative method to use S=3/2 systems as 
a DC vector magnetometer. If the S=3/2 spins fixed 
in a crystal can be rotated around the rotational axis, 
the unambiguous determination of the applied magnetic 
field vector is feasible by monitoring the linewidth of the 
observed ESR spectra while the symmetry axis of the 
crystal is oriented at 9m relative to the rotational axis 
and the rotational axis is moving. This configuration 
also can be realized by producing an array of the crystals 
such that the symmetry axes of each crystal form a cone 
whose opening angle is twice the magic angle. 

These findings provide a better understanding of the 
S=3/2 electronic spin Hamiltonian, especially at low 
fields. They also provide an outlook for the application of 
Vsi in SiC to quantum magnetometry which is promising 
thanks to the electrical properties of SiC, which outstand 
the host material of the NV centers, and the mature fab¬ 
rication technology, which allows an efficient fabrication 
of electronic devices even at the atomic scale^®. 


ACKNOWLEDGMENTS 

We thank Torsten Rendler, Seoyoung Paik, Thomas 
Wolf, Matthias Widmann for helpful discussions, and 
Nathan Chejanovsky for helping to prepare this 
manuscript. We acknowledge funding by the DFG via 
priority programme 1601 and the EU via ERC grant 
SQUTEC and Diadems as well as the Max Planck So¬ 
ciety. 


* s.lee@physik.uni-stuttgart.de 

^ M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, 
J. Wrachtrup, and L. C. L. Hollenberg, Physics Reports 
528, 1 (2013). 

^ R. Schirhagl, K. Chang, M. Loretz, and C. L. Degen, 
Annual Review of Physical Chemistry 65, 83 (2014). 

® E. Sorman, N. T. Son, W. M. Chen, O. Kordina, C. Hallin, 
and E. Janzen, Physical Review B 61, 2613 (2000). 

^ E. Janzen, A. Cali, P. Carlsson, A. Gallstrom, B. Magnus- 
son, and N. T. Son, Physica B: Condensed Matter 404, 
4354 (2009). 

® M. Widmann, S.-Y. Lee, T. Rendler, N. T. Son, H. Fed- 


der, S. Paik, L.-P. Yang, N. Zhao, S. Yang, I. Booker, 
A. Denisenko, M. Jamali, S. A. Momenzadeh, I. Gerhardt, 
T. Ohshima, A. Gali, E. Janzen, and J. Wrachtrup, Nat 
Mater 14, 164 (2015), arXiv: 1407.0180. 

® D. J. Ghristle, A. L. Falk, P. Andrich, P. V. Klimov, J. U. 
Hassan, N. T. Son, E. Janzen, T. Ohshima, and D. D. 
Awschalom, Nat Mater 14, 160 (2015), arXiv:1406.7325. 

^ J. R. Weber, W. F. Koehl, J. B. Varley, A. Janotti, B. B. 
Buckley, C. G. Van de Walle, and D. D. Awschalom, Pro¬ 
ceedings of the National Academy of Sciences 107, 8513 
( 2010 ). 

® W. F. Koehl, B. B. Buckley, F. J. Heremans, G. Calusine, 



and D. D. Awschalom, Nature 479, 84 (2011). 

® P. G. Baranov, A. P. Bundakova, A. A. Soltamova, S. B. 
Orlinskii, I. V. Borovykh, R. Zondervan, R. Verberk, and 

J. Schmidt, Physical Review B 83, 125203 (2011). 

N. Mizuochi, S. Yamasaki, H. Takizawa, N. Morishita, 
T. Ohshima, H. Itoh, and J. Isoya, Physical Review B 
66, 235202 (2002). 

H. Kraus, V. A. Soltamov, D. Riedel, S. Vath, F. Fuchs, 

A. Sperlich, P. G. Baranov, V. Dyakonov, and G. V. As¬ 
takhov, Nat Phys 10, 157 (2014). 

K. Szasz, V. Ivady, I. A. Abrikosov, E. Janzen, M. Bockst- 
edte, and A. Gali, Physical Review B 91, 121201 (2015). 
A. Gali, A. Gallstrom, N. T. Son, and E. Janzen, in Ma¬ 
terials Science Forum, Vol. 645 (Trans Tech Publ, 2010) 
pp. 395-397. 

N. T. Son, P. Carlsson, J. ul Hassan, E. Janzen, 
T. Umeda, J. Isoya, A. Gali, M. Bockstedte, N. Morishita, 
T. Ohshima, and H. Itoh, Physical Review Letters 96, 
55501 (2006). 

R. C. Stevenson, Journal of Magnetic Resonance (1969) 
57, 24 (1984). 

N. M. Atherton, Ellis Horwood series in physical chemistry 
(Ellis Horwood, Ghichester, 1993). 

G. Balasubramanian, I. Y. Chan, R. Kolesov, M. Al- 
Hmoud, J. Tisler, G. Shin, C. Kim, A. Wojcik, P. R. Hem- 
mer, A. Krueger, T. Hanke, A. Leitenstorfer, R. Brats- 
chitsch, F. Jelezko, and J. Wrachtrup, Nature 455, 648 
(2008). 

S. Steinert, Rev. Sci. Instrum. 81, 43705 (2010). 

G. L. Degen, Applied Physics Letters 92, 243111 (2008). 
J. M. Taylor, P. Gappellaro, L. Childress, L. Jiang, D. Bud- 
ker, P. R. Hemmer, A. Yacoby, R. Walsworth, and M. D. 
Lukin, Nat Phys 4, 810 (2008). 

H. Clevenson, M. E. Trusheim, C. Teale, T. Schroder, 
D. Braje, and D. Englund, Nat Phys 11, 393 (2015). 

G. Balasubramanian, P. Neumann, D. Twitchen, 

M. Markham, R. Kolesov, N. Mizuochi, J. Isoya, J. Achard, 
J. Beck, J. Tissler, V. Jacques, P. R. Hemmer, F. Jelezko, 
and J. Wrachtrup, Nat Mater 8, 383 (2009). 

T. Wolf, P. Neumann, J. Isoya, and J. Wrachtrup, ArXiv 
e-prints (2014), arXiv: 1411.6553 [quant-ph]. 

N. Mizuochi, S. Yamasaki, H. Takizawa, N. Morishita, 
T. Ohshima, H. Itoh, and J. Isoya, Physical Review B 
68, 165206 (2003). 

J. Isoya, T. Umeda, N. Mizuochi, N. T. Son, E. Janzen, 
and T. Ohshima, physica status solidi (b) 245, 1298 (2008). 
N. Mizuochi, S. Yamasaki, H. Takizawa, N. Morishita, 
T. Ohshima, H. Itoh, T. Umeda, and J. Isoya, Physical 
Review B 72, 235208 (2005). 

V. A. Soltamov, A. A. Soltamova, P. G. Baranov, and 

I. 1. Proskuryakov, Physical Review Letters 108, 226402 

( 2012 ). 

D. Simin, F. Fuchs, H. Kraus, A. Sperlich, P. G. Bara¬ 
nov, G. V. Astakhov, and V. Dyakonov, Physical Review 


Applied 4, 14009 (2015), arXiv:1505.00176 [cond-mat.mtrl- 
sci]. 

H. Kraus, V. A. Soltamov, F. Fuchs, D. Simin, A. Sperlich, 
P. G. Baranov, G. V. Astakhov, and V. Dyakonov, Sci. 
Rep. 4 (2014). 

J. J. L. Morton, A. M. Tyryshkin, A. Ardavan, 

K. Porfyrakis, S. A. Lyon, and G. A. D. Briggs, The Jour¬ 
nal of Chemical Physics 122, (2005). 

C. Knapp, N. Weiden, H. Kass, K.-P. Dinse, B. Pietzak, 

M. Waiblinger, and A. Weidinger, Molecular Physics 95, 
999 (1998). 

W. Harneit, Physical Review A 65, 32322 (2002). 

S. C. Benjamin, A. Ardavan, G. A. D. Briggs, D. A. Britz, 

D. Gunlycke, J. Jefferson, M. A. G. Jones, D. F. Leigh, 

B. W. Lovett, and A. N. Khlobystov, Journal of Physics; 
Gondensed Matter 18, S867 (2006). 

N. Mizuochi, Y. Ohba, and S. Yamauchi, The Journal of 
Ghemical Physics 111 (1999). 

Y. Teki, S. Miyamoto, M. Nakatsuji, and Y. Miura, Jour¬ 
nal of the American Chemical Society 123, 294 (2001). 

G. Kothe, S. S. Kim, and S. 1. Weissman, Chemical 
Physics Letters 71, 445 (1980). 

Y. Teki, H. Tamekuni, K. Haruta, J. Takeuchi, and 
Y. Miura, Journal of Materials Chemistry 18, 381 (2008). 
J. Isoya, H. Kanda, J. R. Norris, J. Tang, and M. K. 
Bowman, Physical Review B 41, 3905 (1990). 

P. A. van Leeuwen, R. Vreeker, and M. Glasbeek, Physical 
Review B 34, 3483 (1986). 

M. S. de Groot and J. H. van der Waals, Molecular Physics 
3, 190 (1960). 

X. -F. He, N. B. Manson, and P. T. H. Fisk, Physical 
Review B 47, 8809 (1993). 

F. Dolde, H. Fedder, M. W. Doherty, T. Nobauer, 
F. Rempp, G. Balasubramanian, T. Wolf, F. Reinhard, 

L. C. L. Hollenberg, F. Jelezko, and J. Wrachtrup, Nat 
Phys 7, 459 (2011). 

A. L. Falk, P. V. Klimov, B. B. Buckley, V. Ivady, 1. A. 
Abrikosov, G. Calusine, W. F. Koehl, A. Gali, and D. D. 
Awschalom, Physical Review Letters 112, 187601 (2014). 
A. Ajoy and P. Gappellaro, Physical Review A 86, 62104 
( 2012 ). 

E. Bourgeois, A. Jarmola, M. Gulka, J. Hruby, D. Budker, 
and M. Nesladek, ArXiv e-prints (2015), arXiv:1502.07551 
[cond-mat .mes-hall]. 

C. J. Cochrane, P. M. Lenahan, and A. J. Lelis, Applied 
Physics Letters 100, 23503 (2012). 

S. -Y. Lee, S. Paik, D. R. McCamey, and C. Boehme, Phys¬ 
ical Review B 86, 115204 (2012). 

A. Lohrmann, N. Iwamoto, Z. Bodrog, S. Castelletto, 

T. Ohshima, T. J. Karle, A. Gali, S. Prawer, J. C. Mc- 
Gallum, and B. C. Johnson, Nat Commun 6 (2015), 
arXiv: 1503.07566 [cond-mat.mtrl-sci]. 



